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ABSTRACT 

We present a detailed spectral analysis of point-like X-ray sources in the XMM- COSMOS field. 
Our sample of 135 sources only includes those that have more than 100 net counts in the 0.3-10 keV 
energy band and have been identified through optical spectroscopy. The majority of the sources are 
well described by a simple power-law model with either no absorption (76%) or a significant intrinsic, 
absorbing column (20%). The remaining ~ 4% of the sources require a more complex modeling by 
incorporating additional components to the power-law. For sources with more than 180 net counts 
(bright sample), we allowed both the photon spectral index T and the equivalent hydrogen column Nh 
to be free parameters. For fainter sources, we fix T to the average value and allow Nh to vary. The mean 
spectral index of the 82 sources in the bright sample is < V >= 2.06±0.08, with an intrinsic dispersion 
of ~ 0.24. Each of these sources have fractional errors on the value of T below 20%. As expected, the 
distribution of intrinsic absorbing column densities is markedly different between AGN with or without 
broad optical emission lines. We find within our sample four Type-2 QSOs candidates {Lx > 10 44 erg 
s _1 , N# > 10 22 cm -2 ), with a spectral energy distribution well reproduced by a composite Seyfert-2 
spectrum, that demonstrates the strength of the wide field XMM/COSMOS survey to detect these 
rare and underrepresented sources. In addition, we have identified a Compton-thick (Njj > 1.5 x 10 24 
cm~ 2 ) AGN at z=0.1248. Its X-ray spectrum is well fitted by a pure reflection model and a significant 
Fe Ka line at rest-frame energy of 6.4 keV. 

Subject headings: Surveys - Galaxies: active — X-rays: galaxies - X-rays: general - X-rays: diffuse 
background 



1. INTRODUCTION 



Deep penc i l-beam 
(jHasinger et all Il998f ). 



surveys with ROSAT 
Chandra (jBrandt et all 120011 : 



Rosati et all l2002t ICowie et all 120021: lAlexander et all 



20031) and XMM-Newton (jHasinger et aU 120011 : 



Loaring et alJ l2005h have proved that the majority 
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of the X-ray background (XRB) is generated by Active 
Galactic Nuclei (AGN) both in the soft (0.5-2 keV) and 
hard (2-10 keV) band. At fluxes below ~ 10" 14 erg 
cm -2 s _1 in the hard band, the X-ray source population 
in these surveys is mainly comp osed of obscured AGN . 
This supports the suggestion by ISetti fe Woltier I (|1989h 
that the spectral shape of the XRB is due to the inte- 
grated contribution of AGN affected by photoelectric 
obscuration with a wide range of gas column density 
(Nh) and redshifts. Since the resolved fraction of the 
XRB drops fro m » 80 - 90% at 2 - 6 keV to 50 - 70% 
at 6 — 10 keV ([Worslev et alJ |2004[ ) . a sizable number 
of strongly absorbed AGN may still be missing in the 
X-ray surveys. An alternative method to detect heavily 
absorbed AGN is to select objects that have mid-IR and 
radio emission typical of AGN though faint ne ar-IR and 
optical fluxes (Martinez- Sansigre et aTl I2005D . While 
this kind of study cannot quantify which fraction of 
these mid-IR selected, absorbed AGN wou ld be detected 
by X -ray selection, the COSMOS survey (|Scoville et al.1 
l2007h will be able to answer this question due to its rich 
multi-wavelength coverage (from radio to X-ray) on a 
large area of the sky (2 deg 2 ). 

The XMM-Newton wide- field survey in t he C OSMOS 
field (XMM-COSMOS, lHasinger et aU 127)071 ). with 
an unprecedented combination of wide area coverage 
and high sensitivity, is providing a large number of 



2 



Mainieri et al. 




Fig. 1. — Net [0.3-10] keV pn counts distribution for the sample 
of 135 X-ray sources used in this work. 

AGN with enough counts to perform a detailed study 
of their X-ray spectra. This spectral information, 
particularly the Nh distribution, is a fundamental input 
parameter to mo del the XRB (e.g. IComastri et al.lll995t 
iGilli et aLll2001l ). While we anticipate the completion of 
the multi-wavelength campaigns including the optical 
spectroscopic follow-up within the next few years, we 
report in this paper the X-ray spectral fitting results 
for a preliminary sample of spectroscopically-identificd 
X-ray sources. The paper is structured as follows: in §2 
we describe the sample selected on the basis of counts 
statistics and optical identification; in §3 we describe our 
X-ray spectral extraction procedure, in §4 we present 
the results of the X-ray spectral analysis, in §5 we 
discuss the properties of four Type-2 QSOs, in §6 we 
compare the X-ray and optical classification, and finally 
we summarize our conclusions in §7. 
Throughout the paper we assume Hq = 70 km s _1 
Mpc" 1 , n m =0.3 and Ov=0.7. 

2. SAMPLE SELECTION 

XWLM-Newton has imaged the full 2 deg 2 of the COS- 
MOS area down to the following flux limits in the respec- 
tive energy bands: 7 x 10~ 16 erg cm -2 s" 1 [0.5-2 keV], 
4.0 x 10" 15 erg cm~ 2 s" 1 [2-10 ke V] and 1.0 x lO"* 4 erg 
cm" 2 s- 1 [5-10 keV] (see Fig. 7 of lCappellutilFaT1l2C^ 
for details on the sky coverage as a function of the X-ray 
flux). A general out line of the survey can be found in 
iHasinger et alJ (|2007h . Further details such as the point- 
source detection method and sky area c overage as a func- 
tion o f the X-ray flux are presented in ICappelluti et all 
(|2007h . Our sample is b ased on the X-ray cata logue 
of 1390 point-like sources (|Cappelluti et alJl200l . We 
limit our analysis to the sources detected w ith the EPIC 
pn— CCD (pn) camera (jStruder et al .1 l200lh . in the first 
12 XMM-COSMOS observations sinc e optical sp ectro- 
scopic follow-up (|Trump et all 120071 : iLillv et all l2007h 
has been concentrated in this area (~ 1.3 deg 2 ). Th ese 12 
fields are flagged in Table 1 oflHasineer et al.l (120071) . Re- 
liable optical counterparts (jBrusa et al.ll2007l ) have been 
determined for ~ 90% of the sources in these 12 fields. 



Fig. 2. — X-ray [0.5-10 keV] flux distribution for all the X-ray 
sources (empty histogram) and for the sample of spectroscopically 
identified sources (filled histogram), with more than 100 net counts, 
we analyze in this work. 

We exclude 20 of the 715 X-ray sources in this area that 
are classified as "extended" from the detection algorithm. 
The observed X-ray emission from these sources is likely 
to be due to a group or cluster of galaxies, while here 
we are interested in selecting AGN. From the remaining 
695 X-ray sources, we select sources with greater than 
100 net counts in the [0.3-10] keV energy band and opti- 
cal spectroscopic identification. We furthe r remove one 
sourc e that has been identified as a star (|Trump et al.l 
120071 ). The final sample comprises 135 objects. We show 
the distribution of their net counts in the [0.3-10] keV 
band in Figure [T] and the [0.5-10] keV flux distribution 
(FigurcEJ) that covers a range of 1.4 x 10~ 15 , 1.2 x 10~ 13 
erg cm~ 2 s _1 . From their optical spectra, we can fur- 
ther subdivide our sample based on the presence of broad 
emission lines: 'Broad Line AGN' (BL AGN, 86 objects; 
FWHM>2000 km s" 1 ), 'NON-Broad Line AGN' (NON- 
BLAGN, 49 objects; FWHM<2000 km s" 1 ). We note 
that in this latter class there are objects showing clear 
signs of nuclear activity such as high excitation emission 
lines, as well as sources with normal galaxy spectra. We 
compare this purely optical classification with the X-ray 
properties of our sources in Sec. [6] 

3. EXTRACTION OF X-RAY SPECTRAL 
PRODUCTS 

We have implemented an automated procedure to pro- 
duce the X-ray spectrum for each source by combining 
counts from individual exposures. We have used the lat- 
est release of the XMM- Newton Science Analysis System 
(SAS) 16 software package (v 7.0). The task region has 
been used to generate the source and background ex- 
traction regions. The source region is defined as a cir- 
cle with radius r s that varies according to the signal-to- 
noise and the off-axis angle of the detection to optimize 
the quality of the final spectrum. The radii of these re- 
gions are reduced by the task to avoid overlapping with 
the extraction regions of nearby sources. All source re- 

16 http:/ /xmm. vilspa.esa.es/external/xmm_sw_cal/sas_frame.shtml 



X-ray spectral properties of AGN. 
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Fig. 3. — Le/4: spectral slope value (V) from the fit of each single source using the PL model as a function of the net counts in the full 
[0.3-10] keV band. Right: fractional statistical error (lcr) on T as a function of the net counts in the [0.3-10] keV energy band. The dashed 
line in both plots is the threshold of 180 net counts in the [0.3-10] keV band that divides the bright from the faint sample. 



gions are further excised from the area used for the back- 
ground measurement. The task especget has been used 
to extract from the event file the source and background 
spectra for each object. The same task generates the 
calibration matrices (i.e. arf and rmf) for each spectrum 
and determines the size of the source and background 
areas while updating the keyword BACKSCAL in the 
header of the spectra appropriately 17 . The single point- 
ing spectra have been combined with mathpha to gener- 
ate the spectrum of the whole observation. 18 For each 
source in our sample, we use all the available counts from 
the XMM- COSMOS observations, including those com- 
ing from over lapping fields not inc luded in the 12 fields 
list (see Fig.l lHasinger et aLll2007|) . Finally, in order to 
use the x 2 minimization technique, we bin the spectra 
with grppha so that each bin has at least 20 counts. 

4. SPECTRAL ANALYSIS 

iTozzi et alJ (|2006| ) have shown by extensive simula- 
tions that below 50 co unts the bes t fit values obtained 
using Cash statistics (ICash I Il979t) are more accurate 
than those obtained with the % 2 . For greater than 50 
counts, the two methods give equivalent results. Since 
we limit our analysis in this paper to sources with more 
than 100 counts, we are confident that the results ob- 
tained with the x 2 minimization technique are accu- 
rate. We use XSPEC 19 (vll.3.2p) for our spectral fit- 
ting analysis. We first fit the data with two basic in- 
put models: a simple powerlaw (PL) and a powerlaw 
modified by intrinsic absorption at the redshift of the 
source (APL). Both models include an additional com- 
ponent to account for photoelectric absorption due to 
the Galactic column density that is fixed to the value 
in the COSMOS region (N§ al ~ 2.7 x 10 20 cm- 2 , 
iDickev fc Lockman I [l990h 20 . The PL model is made 

17 The header keyword BACKSCAL is set to 1 for the source 
spectrum while for the background spectrum it is fixed to the ratio 
between the background to source areas. 

18 We note that all the XMM-Newton observations in the COS- 
MOS field have been performed with the thin filter for the pn 
camera. 

19 http:/ /heasarc. gsfc.nasa.gov/docs/xanadu/xspec/ 

20 This is an average value for the Galactic N H in the COSMOS 
area where N§ al is in the range [2.5-2.9] xl0 20 cm -2 . This range in 



of two XSPEC components wabs*zpowerlw, while the 
APL model consists of the combination of three differ- 
ent components wabs*zwabs*zpowerlw. The wabs model 
describes the photoelectric absorption using W isconsin 
cross-sections (|Morrison fc McCammonl 11983ft and its 
only parameter is the equivalent hydrogen column den- 
sity (zwabs has the redshift as an additional parameter). 
zpowerlw is a simple power law parameterized by the pho- 
ton index, the redshift and a normalization factor. 21 The 
model fits yield the power-law photon index T, the X- 
ray luminosity in the [0.5-2] and [2-10] keV rest-frame 
bands, and from the APL model also the intrinsic col- 
umn density Nh ■ We notice that the dispersion of L for 
our sample increases significantly for sources with low 
counts statistics (Fig. [3l left panel) and in particular 
the fractional error becomes quite large (Fig. 02 right 
panel). Above 180 net counts, the fractional error re- 
mains below 20%. Hence, we split our sample in two: 
sample-1 including 82 sources with more than 180 net 
counts, and sample-2 having 53 sources with less than 
180 counts. For sample-1, we allow both T and Nh free 
to vary, while we fix T to the average value, obtained 
with sample- 1, for lower count sources (sample-2). 
For all the 135 X-ray sources, we perform a spectral fit 
using both PL and APL models. We label a source as 
X-ray absorbed in those cases for which the APL model 
is a better fit, than the pure PL, with a confidence level 
threshold of 90% based on an F-test. 
The output of our spectral analysis is reported in Tab. 
[H The table has the following structure: IAU name 
(col.l), identification number (xid, col. 2), X-ray coordi- 
nates (col. 3-4), net detected X-ray counts in the [0.3-10] 
keV band (col. 5), spectroscopic redshift (col. 6), best 
fit model (col. 7), spectral index T (col. 8), intrinsic 
column density Nh (col. 9), X-ray fluxes (col. 10-11-12), 
de-absorbed X-ray luminosities (col. 13-14-15). 

4.1. Notes on some individual sources 

Galactic column density does not affect the results of our spectral 
analysis. 

21 We refer the reader to 

http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/ for further 
details on the spectral models. 
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TABLE 1 

Parameters of the best fit model for sources with soft excess 



XID counts a Model b x 2 /d.o.f. T c N H d kT c T f Redshift Opt. class 8 



41 


315 


APL+po 
APL+bb 


0.94 
0.83 


i-95?:ts 


oi SS 21 - 59 
z±.oo 21 .0 2 

91 ^1 21.55 
z±.o± 20 42 


30±§ 


2.0 


0.114 


NLAGN 


106 


141 


APL+po 
APL+bb 


0.25 
0.26 


2.0 
2.0 


99 oo22.66 
"' 00 21.98 
99 9Q.23.00 
ZZ.Z0 21 gg 


i2ii 7 j 4 


2.0 


0.710 


gal 


117 


111 


APL+po 


0.69 


2.0 


22.76| 2 ;3| 


01 +48 
Si -18 


2.0 


0.936 


gal 






APL+bb 


0.59 


2.0 


99 c;v22.86 
zz - o, 20.42 








274 


112 


APL+po 
APL+bb 


0.26 
0.24 


2.0 
2.0 


99 fi'723.00 
ZZ.D( 22 lg 

99 CC23.00 

zz.oo 2 q 42 


1 24+ 878 
1Z4 -119 


2.0 


0.677 


gal 



a Net pn counts in the [0.3-10] keV energy range. Best fit model: APL+po— absorbed power-law plus an extra power-law 
for the soft excess; APL+bb— absorbed power-law plus a blackbody for the soft excess. c Slope of the powerlaw model 
(photon index). Logarithm of the intrinsic absorption (cm _2 ).° Temperature (in eV) of the blackbody used to model 
the soft excess/ Slope of the extra power-law used to model the soft excess. s Optical classification, see Sec. [6] for details. 



For each source we carefully check the results obtained 
with the basic PL and APL models and, if significant 
residuals are present, we refine the fit using more com- 
plex models. We show in Fig. [4] a representative X-ray 
spectrum for each one of the different best-fit models. We 
use the F-test and a confidence level threshold of 90% to 
choose between the different models. 

4.1.1. Soft excess 

A clear soft excess is present in four of our sources 
(xid 41, 10 6, 117, 274). This feature, first observed wit h 
EXOSAT (jArnaud et al.lll985l : iTurner fc PoundslflQSQl ). 
has been confirmed by XMM-Newton observations (e.g. 
Pounds fe Reeves 1120021 : iPorquet et aHl2004t iGaTlo et all 
200(f ) . but its origin is still uncertain. Such a soft com- 
ponent may be the high energy tail of the UV bump 
(a blackbody model is appropriate in this case), or can 
be due to reprocessed emission scattered along our line 
of sight by a photo-ionized gas located just above the 
obscuring torus (an additional power-law with the spec- 
tral index fixed to the value of the hard X-ray primary 
power-law is a good parameterization of this scenario). 
We fit these four sources adding to the basic APL model 
an extra component represented either by a power-law 
(po) or a blackbody (bb 22 ) according to the two physical 
scenarios mentioned above. We report in Tab. [T]the pa- 
rameters of the additional component in the fit of these 
four sources. We are not able to distinguish on a sta- 
tistical basis between the two models given the similar 
values of x 2 . Nevertheless, we notice that all these four 
sources present intrinsic absorption and therefore we ex- 
clude that the soft-excess of these objects is due to the 
high energy tail of the UV bump (APL+bb). 

4.1.2. Fe Ka line 

Three sources show significant features ascribable to 
the rcdshifted Fe Ka emission line: xid 2028, 2043, 2608. 
For these sources we add a Gaussian component (gauss) 
to the model, fixing the line energy to 6.4/(1 + z) keV. 
The best fit values of interesting parameters are reported 
in Tab. [H We show in Fig. [5] the ratio of the data versus 
the model (powerlaw for xid 2028, 2043 and pexrav 23 for 

22 This is a blackbody spectrum defined by the temperature kT 
in keV and a normalization factor. 

23 An exponentially cut off power law spectrum reflected from 
neutral material. We refer the reader to Magdziarz & Zdziarski 



TABLE 2 

Parameters of the gauss additional component 

FOR THE SOURCES WITH Fe LINE 



XID 




EW b 


Redshift 


Opt. class 


2028 
2043 
2608 


61 6+ 364 
u±u -224 
m +120 
i,y -115 
901 +408 


2754ti«» 
748^ 
792t 4 gf 


0.784 
0.668 
0.125 


gal 
gal 
gal 



a Observed width of the line in eV. Rest frame equivalent 
width of the line in cV c Optical classification, see Sec. [6] for 
details. 



xid 2608) in an energy range around the expected loca- 
tion of the Fe Ka line. Interestingly all three of these 
sources do not show sign of AGN activity from their op- 
tical spectra and are therefore classified as 'galaxy'. 

4.1.3. Thermal emission? 

Source xid 54, if fitted with an APL model, gives a 
large value for the spectral slope (r > 3) and signifi- 
cant residuals in the 0.3-10 keV energy range. An alter- 
native description of its spectrum is obtained assuming 
we are observing therma l emission, parameterized w ith a 
Raymond-Smith model (|Ravmond fc Smith Ill977[) with 
a temperature kT= 1.6ig 2 keV fixing the metalicity to 
0.3 solar. Source xid 54 is identified with two interacting 
galaxies (see Fig. [6|) at redshift z= 0.350 with no sign 
of AGN activity from its optical spectrum. Its X-ray lu- 
minosity of 3 x 10 42 er g s -1 is larger than tha t expected 
for early-type galaxies (jMatsushita et al.l2001l ) and, from 
the optical imaging, there is a concentration of galaxies 
around xid 54 with the same photometric rcdshifts there- 
fore supporting the idea that we are looking at the X-ray 
emission from a group of galaxies. Nevertheless, we can 
not exclude with the current data that a fraction of the 
X-ray flux of source 5 4 could come from an absorbed nu- 
cleus (e.g XBONGS. IComastri et al1l200l or from dis- 
crete sources like LMXBs or HMXBs in the galaxy. A 
Chandra observation with its higher angular resolution 
could possibly locate discrete sources inside xid 54. 

4.1.4. XID=2608 : a Compton-thick AGN? 
(1995) for a detailed description of such model. 



X-ray spectral properties of AGN. 



XID= 1 z = 0.373 



XID=293 z = 0.444 




channel energy (keV) 



channel energy (keV 



Fig. 4. — Examples of X-ray spectra with different best fit model. Top-left: unabsorbed power-law (PL); top-right: absorbed power- 
law (APL); bottom-left: absorbed power-law plus a black-body component to model the soft excess; bottom-right: thermal emission 
parameterized with a Raymond-Smith model. 



An additional source that requires a more complex 
modeling of its spectrum is xid 2608. The fit with the 
APL model gives an extremely flat value for T (« 0.3) 
and large residual s at both low an d high energies (see left 
panel of Fig. [7}. iHasinger et all (|2007l ) found that this 
source is located in an area populated by local Compton- 
thick Seyfert-2 galaxies in an X-ray co lor-color diagram 
(see Fig. 12 of lGuainazzi et al.ll2005l) . This, together 
with other evidences based on lines ratios from the opti- 
cal spectrum, supports the hypothesis that source 2608 
is a heavily absorbed AGN. 

We use the 131 net counts from the pn camera for this 
source to study more in detail its X-ray spectrum. A 
pure reflection component model (pexrav) is a better de- 
scription compared to the APL model according to an 
F-test with a confidence level of 95%. Nevertheless this 
fit leaves a clear residual around the expected position of 
the 6.4 keV Fc Ka line. The best-fit model for xid 2608 
is a pure reflection model plus a Gaussian line at 6.4 keV 
rest- frame (pexrav + gauss). The details for the differ- 
ent spectral fits are reported in Tab. [3J The presence 
of the Fe Ka fluorescent line at 6.4 keV is significant at 
95% according to an F-test. The presence of the line is 
a clear sign that the source is heavily absorbed, but a 
useful observable to confirm its Compton-thick nature is 
the Equivalent Width (EW) of the same line. The nom- 
inal best fit value for the EW (792T^gf eV) is higher 
than the maximum (600 eV) observed EW in Compton- 



thin objects (|Turner et al.| [l997). This supports the idea 
that source 2608 is a Compton-thick AGN, although we 
have to mention that with the current photons statistics, 
the error for the observed flux of the line ( and conse- 
quently for the EW) is still large. We are confident that 
an improved result will come after the completion of the 
additional 600 ksec XMM observations awarded in A04. 
Another diagnostic on the Compton-thick nature of 
this source could be the thickness parameter T=F[2- 
10keV]/F[OIII]. A high quality optical spectrum for this 
source is available in the Sloan Digital Sky Survey archive 
and we obtain a value fo r F[OIII] from the analysis of 
iKauffmann et al. I (|2003h . The [OIII] flux has bee n cor- 
rected for the extinction toward the narrow-line region as 
deduced from the Balmer decrement. We obtain T= 3.8 
which is in a "grey area" where both Compton thick and 
le ss absorbed AGN ar e located (see for example Fig. 1 
of lBassani et al.lll999D . 

In Fig. 1131 we assume as a lower limit for the column 
density of xid 2608 the value 1.5 x 10 24 cm -2 where the 
Compton optical depth is equal to unity and the directly 
transmitted nuclear emission is strongly suppressed in 
the [0.3-10] keV band. For the luminosity of this object, if 
we assume that only 3% of the flux has been reflected, we 
obtain a value of ~ 7.4 x 10 43 erg s -1 , while for reflected 
fractions between 10— 1% the luminosity would be in the 
range (0.2 



2.2) x 10 44 erg s" 1 . 
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XID=2028 z=0.784 



TABLE 3 

Parameters of the best fit model for source xid 
2608 



Modcl a 



EW C 



d.o.f. 



channel energy (keV) 
data/ mode 



XID=2043 z=0.668 



channel energy (keV) 



XID=2608 z=0.125 



APL 

pexrav 
pexrav+gauss 



2.0 
2.0 
2.0 



0.16: 



792Z 



9.3 
4.1 
1.7 



11 
9 



a Best fit model: APL — absorbed power-law; pexrav — pure 
reflection model; pexrav+gauss — pure reflection model plus 
a Gaussian line. Hydrogen column density in unit of 10 22 
cm -2 . c Equivalent width of the Fe Kq line expressed in eV. 




Fig. 6. — ACS image of source xid 54. The cutout is 90 arcsec 
on a side. 



channel energy (keV) 

Fig. 5. — The ratio of the data versus a powerlaw model (top 
panels) or a pexrav model (bottom figure) around the energy of 
the Fe Ka line for the three sources with significant detection of 
this feature. 



4.2. Spectral properties of the sample 

As mentioned in Sec. [4] we leave both T and Nh free 
to vary when fitting the sources in sample-1. The results 
of this analysis are summarized in Fig. [5] The average 
value of r docs not change as a function of Nh as a lready 
noticed in deep surveys (i.e. iMainieri et al.l l2002f). We 
obtain, using the weighted mean, < T >= 2.06 ± 0.08 
and the observed dispersion of the distribution of the 



best fit values is a 0.25. As the typical error in a 
single measurement of T is Ar = 0.09, assuming that 
both statistical errors and the intrinsic dispersion are 
distributed as a Gaussian, the intrinsic scatter in T is 
Cint ~ 0.24. For comparison with X-ray sp ectral studies 
in a sim ilar X-ray flux range of our sample, iMateos et al.l 
(|2005a|) from a large sample of serendipitous sources de- 
tected with XMM-Newton in a ~ 3.5 deg 2 area, ob- 



tained < T >= 1.96 ± 0.01; iPerola et all (|2004h in the 
spectroscopic analysis of the HELLAS2 XMM ldf found 
< T >= 1.90 ± 0.22; iPage et all (|2006f ) from the spec- 
tral fit of AGN in the 13 H XMM-Newton/ Chandra deep 
field found < T >= 2.0 ±0.1 with an intrinsic dispersion 
a w 0.36. All these measurements are consistent with 
each other within the uncertainties. Adopting the optical 
classification described in Sec. [2] the mean value for the 
spectral slope for BLAGN (58 sources) is < T >= 2.09 
with a dispersion of a ss 0.26, while for not BLAGN (24 
sources) we obtain < V >= 1.93 and a w 0.29. Further- 
more, we confirm that the average value of the photon 
index does not vary with redshift in the range z= [0.0,3.0] 
covered by our sa mple, thus confirming p revious findings 
(e.g. see Fig. 9 of iPiconcelli et all (|2003f ) for a compila- 
tion from the literature) 

The other physical quantity that we measure from the 
spectral fitting is the column density Nh- In this case, 
we consider all our 135 sources since Nh has been left 
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channel energy (keV) channel energy (keV) 

Fig. 7. — The X-ray fit of source xid 2608 with the basic APL model (left) and a pure reflection model plus a Gaussian line (right). 



free to vary both in sample-1 and sample-2. According 
to an F-tcst, 32 X-ray sources do require intrinsic ab- 
sorption in excess to the Galactic one, at a confidence 
level larger than 90%. Therefore a fraction as large as 
24% of our sample is made of X-ray absorbed AGN. Fig. 
[9] shows the distribution of Nh values for these sources. 
We note that the observed Nh distribution refers only 
to the sources inside the region in the Nh-Lx-z space 
delimited by the count-rate detection threshold of our 
survey. This introduces a bias against absorbed sources, 
and therefore the fraction of absorbed sources detected 
in our study has to be considered a lower limit. Ac- 
cording to the most re cent p opulation synthesis model 
of the XRB (jGilli et al.ll2006f ). in the band used to select 
our sample the expected fraction of obscured source with 
column density Nh > 10 21 cm -2 is 20% that is consis- 
tent with what we found. In Fig. [9] we divide BL AGN 
from NOT BL AGN. The visual impression that NOT BL 
AGN have larger column densities than BL AGN is con- 
firmed by a Kolmogorov-Smirnov test that gives a prob- 
ability larger than 99.9% that the two distributions are 
different. No object with Nh > 10 22 cm -2 shows broad 
lines in its optical spectrum. Nevertheless, 9% (8/86) 
of the BL AGN in our sample do show s ome intrinsic 
absorption in their X-ray spectra (see al s o iMittaz et aT 
1999 HFiore et al.| l200l l: [Page et al.l |200H ISchartel et al 
Tozzi et alj|2001t [Mainieri et alj |2002t i Brusa et al 



2001] 

200a 



Perola et aij|2004 iMateos et alJl2005blT 



In Fig [10] we show the distributions of the R-K (Vega) 
colors for the sources with 'PL' as best fit model (empty 
histogram) and for the ones that instead require an 
absorbed power-law (hatched histogram). The two 
distributions are significantly different according to a 
Kolmogorov-Smirnov test with a probability of 99.99%. 
The X-ray sources that require an absorption compo- 
nent in their spectral fit are on average redder, suggest- 
ing a correlation between X-ray absorption and optical 
to near-IR colors. On the contrary, sources that do not 
show absorption in their X-ray spectra have bluer color 
typical of optically selected, unobscured quasars. These 
results confirm those obtained from an a nalysis based on 
HR values made by IBrusa et alJ ()2007l ) (see their Fig. 
10). Nevertheless the interpretation of this correlation 
between X-ray absorption and optical to near-IR colors is 
not straightforward since we are sampling different scales 



in the two measurements (i.e. nucleus with the X-ray 
data and nucleus+host galaxy with the R-K colors). 

When the number of counts in a source is inadequate 
to perform a spectral fit, a widely used tool to study 
the general spectral properties of an X-ray source is the 
hardness ratio HR = (H — S)/ (H + S), where H are the 
counts in the [2-4.5] keV band and S those in the [0.5-2] 
keV energy band. In Fig. [TTJwc show the hardness ratio 
values versus the amount of intrinsic absorption derived 
from our spectral analysis (both sample-1 and sample- 
2). A clear correlation between the two quantities is 
present: 90% 24 of the sources with Nh > 10 22 cm~ 2 
have HR> -0.3 and 99% of the sources with N H < 10 22 
cm~ 2 have HR< —0.3. Therefore, although one has to 
remember tha t the HR is a stron g function of redshift 
(e.g. Fig. 8 in lSzokolv et al.|[2004l ). it is still possible to 
use HR for statistical studies. 

Another diagnostic that can yield important infor- 
mation on the nature of X-r ay sources is the X- 
ray-to-optical flux ratio (e.g., iMaccacaro et al.1 119881 ; 
iStocke et ail |1991[) . The majority of the AGN have 
X-ray -to-optical flux ratios (X/O) of 0.1 <X/Q < 10 
(e.g., lAkivama et~afl l2000t iLehmann et~afl l200lt ). but 
Chandra and XMM-Newton surveys have shown that 
there is a non negligible population of AGN with high 
X/O (> 10) and that a large fraction of them are ob- 
scured, and p ossibly high-redshift, Type-2 QSOs (e.g. , 



. p 

Fabian et al. 2000; Mainieri et al. 2002; Fiore et al.H2003l : 



Mignoli et "all 120041 : iMainieri et"aTll2005l ). For compari- 
son with the literature, we define X/O as the ratio be- 
tween the X-ray flux in the [2-10] keV band and the flux 
in the optical R band. In Fig. [12] we plot the X/O values 
for the sources in our sample versus the Nh- Out of the 
seven sources that have X/0> 10, four show absorption 
in their X-ray spectra (APL) and one is a Typc-2 QSO. 
We notice that the other three Type-2 QSO candidates in 
our sample (see Sec. [5]) have X/O values inside the range 
0.1 <X/0< 10 where most of the optical or soft X-ray 
selected AGN are located. Since we limit our analysis to 
the brighter X-ray sources and the spectroscopic follow- 
up is not complete, we postpone any further analysis on 
the nature of X/0> 10 sources to a future paper. 



We note that the only source with Nh > 10 
HR< —0.3 shows a soft excess in its X-ray spectrum. 
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Fig. 8. — r versus Nh for the X-ray sources with more than 180 
net counts in [0.3-10] keV (bright sample) and spectroscopically 
identified. Filled circles are BL AGN, while empty circles are not 
BL AGN. Error bars correspond to lcr. To simplify the figure, 
we did not report the error bars on T for unabsorbed sources and 



plotted them to N H =N 



gal 



2.7 X 10 20 cm 




21 22 23 
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Fig. 9. — Intrinsic column density (Nh) distribution for BL AGN 
(empty histogram) and NOT BL AGN (hatched histogram) with 
intrinsic absorption in excess of the Galactic column density. 



Fig. 10. — R-K colors (Vega) distribution for sources with PL as 
best fit model (empty histogram) and for sources with 'APL' as a 
best fit model (hatched histogram). 
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Fig. 11.— HR defined using the [0.5-2] (soft) and [2-10] (hard) 
bands versus the column density derived from the spectral fitting 
analysis. Only sources with errors on the HR smaller than 0.3 have 
been plotted. Filled circles are BL AGN, while empty circles are 
NOT BL AGN. The horizontal dashed line corresponds to HR= 
—0.3 used to separate absorbed and unabsorbed sources, while the 
vertical dashed line indicates a column density equal to 10 22 cm -2 . 



5. TYPE-2 QSO CANDIDATES 

Using the spectral parameters from the best-fit model, 
we correct the X-ray luminosity of each source for the 
intrinsic and Galactic absorption. These corrected lumi- 
nosities are plotted in Fig. [13] versus the Nh for all the 
sources in our sample. Four objects are characterized by 
a high X-ray luminosity (Lx[0.5 — 10 keV] > 10 44 erg 
s" 1 ) and substantial absorption (Nh > 10 22 cm -2 ) and 
we can therefore classify them as Type-2 QSOs. Radio- 
loud Type-2 QSOs a re known since l ong times thanks 
to radio surveys (see iMcCarthvl 119931 for a comprehen- 
sive review), while radio-quiet Type-2 QSOs have been 
observed only recently in Chandra and XMM- Newton 
X-ray surveys (|Dawson et "al] [200ll ; iNorman et all [20021 ; 



Alain 



eri et all [2002; Stern et al]l200l iDella Ceca et all 



20031 iFiorc et al. 2003; Tozzi et al|2006ft and optical sur- 
veys (SPSS. IZakamska et al J 1 2003 ). Two of our Type-2 
QSOs candidates, xid= 70,2289, are clearly detected in 
the radio at 20 cm using the Very Large Array (VLA) 
with an integ rated flux of 540 ± 24 a nd 52 ± 11 microJy 
respectively (jSchinnerer et all |2007[ ) . The radio power 
of these two sources is therefore Pi.4GHz = 9.8 x 10 23 
and 1.5 x 10 23 W/Hz. Historically such radio power has 
been used to divide radio loud and radio quiet AGN 
but such a dividing line appears to be rcdshift depen- 
dent: « 5 x 10 23 W/Hz for the Palomar Green sample 
(mainly below z < 0.3) up to 5 x 10 25 W/Hz for the 
Large Bright Quasar Survey sample (< z >^ 1.2). Since 
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Fig. 12.— X/O (f[2-10 keV]/F(R)) versus N H values. The filled 
circles are BL AGN, the empty circles NOT BL AGN. We label the 
four Type-2 QSOs candidates. The horizontal dashed line indicates 
the value X/O=10. 

our two sources are at z~ 0.7 — 0.8, we suggest to clas- 
sify them as radio quiet AGN. The other two objects 
(xid= 122, 2237) are not detected in the radio and we 
can fix a 4. 5 a upper limit to their radio flux of 50 and 
54 micro Jy (|Schinnerer et alj|20"07t ). 
The optical spectra of these four sources show high exci- 
tation emission lines and their redshifts are 0.688, 0.831, 
0.941, 2.418 respectively for xid 70, 2289, 2237, 122. 
Using the multi-b and photometry a vailable from the 
COSMOS survey (jCapak et al.ll2007h . we have derived 
the spectral energy distribution (SED) for the four 
Type-2 QSOs and compared them with the spectrum of 
NGC6240 and a Seyfert-2 compo site spectrum derived 
fro m a sample of l ocal g alaxies by iSchmitt et al.l (|1997f ) 
and iMoran et all (pOOll ). While the SED of NGC6240 
does not reproduce well the observed photometry of our 
Type-2 QSOs, an excellent description of the same is 
given by the composite Seyfert-2 SED (see right panels 
in Fig. [14)) . Furthermore, the R-K colors of these four 
objects are red (R-K=4.58, 3.91, 4.97, 4.76 respectively) 
although they can not be classified as EROs (R-K> 5). 

6. COMPARISON BETWEEN X-RAY AND 
OPTICAL CLASSIFICATIONS. 

A classification based on the properties of the optical 
spectra of the 135 sources in our sample divides them 
into 'Broad Line AGN' (BLAGN, 86 objects) if emis- 
sion lines broader than 2000 km s _1 are present, 'Narrow 
Line AGN' (NLAGN, 32 objects) if the optical spectrum 
shows high excitation emission lines and 'galaxy' (gal, 
17 objects) if there is no sign of AGN activity from the 
optical spectrum. As shown by deep Chan dra and XMM- 
Newton surveys (e.g. ISzokolv et al.ll200"4[ ) a pure optical 
classification of AGN is biased against absorbed sources 
that appear as nor mal galaxies at those wav elengths. As 
previo usly done by ISzokolv et al.1 (|2004l ) and iTozzi et ail 
pool , we introduce an X-ray based classification: we 
define X-ray absorbed AGN sources that are best fitted 
by an APL model compared to the PL one and have 
Lx > 10 42 erg s _1 , X-ray unabsorbed AGN sources best 



Fig. 13. — Intrinsic, de-absorbed X-ray luminosity in the [0.5- 
10] keV band vs. N H . The filled symbols are BLAGN, while the 
empty symbols are not BLAGN. For source xid 2608 we assume a 
lower limit on Njj of 1.5 X 10 24 cm -2 and for the luminosity we 
estimate a value of ~ 7.4 X 10 43 erg s -1 assuming that a fraction of 
3% is reflected (the error bar shows the luminosity range covered 
assu ming t hat the reflected fraction is between 1% and 10%). See 
Sec. 14.1.41 for details. The dashed lines define the "Type-2 QSO 
region" . 

fitted with a PL model and Lx > 10 42 erg s _1 and finally 
X-ray galaxies sources with Lx < 10 erg s^ 1 . 
Table [5] shows the comparison of the optical and X-ray 
classifications for our 135 sources. Ninety-one of these 
sources (~ 67%) have a similar classification from the 
optical and X-ray data. The best agreement between 
the two classifications is for Broad Line AGN (optical) 
and X-ray unabsorbed AGN (X-ray) for which the frac- 
tions of similar classifications are of the order of 91% 
(78/86 Broad Line AGN) and 76% (78/102 X-ray unab- 
sorbed AGN) respectively. The - 9% of BL AGN that 
show X-ray absorption in their X-ray spectra have values 
of the column density Nh below 10 22 cm~ 2 (see empty 
histogram in Fig. [9j. The main difference is instead for 
objects classified as galaxies on the basis of the optical 
spectra. Most of these objects (16/17) are classified as 
AGN (11 absorbed and 5 unabsorbed) on the basis of the 
X-ray luminosity. This confirms that the X-ray classifi- 
cation is more successful than the optical one in revealing 
the presence of black hole activity. The situation is inter- 
mediate for Narrow Line and X-ray absorbed AGN: only 
~ 41% of the optically classified Narrow Line AGN do 
show detectable X-ray absorption. We note that of the 
remaining Narrow Line AGN, ~ 80% have z> 0.4 and 
therefore the Ha line is outside the observed wavelength 
range, while for nine of them the Mgll line is inside the 
observed range (i.e. 0.92 < z < 2.29) but the S/N of the 
spectra could not be sufficient to detect a weak broad 
line. It is therefore possible that at least part of the dis- 
agreement between the optical and the X-ray classifica- 
tions for these objects is due to less than optimal optical 
spectra, in terms of either spectral coverage or S /N. 

7. CONCLUSIONS 

We have presented the detailed spectral analysis of 135 
X-ray sources from the X.MM.- Newton wide-field survey 
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Fig. 14. — Left: i-band (F775W) ACS cutouts for the four Type-2 QSO candidates. Each cutout is 10 arcsec across. Right: the spectral 
energy distribution of the Type-2 QSO candidates (filled circles) compared to the SED of a composite Seyfert-2 spectrum (dashed line) 
and NGC6240 (dotted line). 



in the COSMOS field. All the sources in our sample 
have more than 100 net counts in the [0.3-10] keV en- 
ergy band and have been spectroscopically identified. For 
each source we have performed an accurate spectral fit 
in order to measure the continuum shape, the amount 
of absorbing matter and the strength of other spectral 
features. Our main results are summarized as follows: 

• We find that, to the X-ray flux limit we are sam- 



pling (F x [0.5 - 10] = 1.4 x 10" 15 erg cm- 2 s" 1 ), 
~ 76% of the spectra are well reproduced with a 
single power-law model, ~ 20% require an absorbed 
power-law model and the remaining ~ 4% need 
more complex models. 

• The average value of the spectral slope of the in- 
trinsic spectrum for the 82 sources with more than 
180 net counts {Sample-!) is < Y >= 2.06 ± 0.08 
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with an intrinsic dispersion of cr,. 



0.24. 



• We find no correlation between the spectral slope 
r and the amount of intrinsic absorption Njj, con- 
firming that the hardening of the X-ray spectra go- 
ing to fainter X-ray fluxes is due to the increased 
fraction of absorbed X-ray sources. 

• None of the X-ray sources with a column density 
Nh > 10 22 cm~ 2 shows broad line in their opti- 
cal spectra, although a fraction (9%) of broad line 
AGN shows intrinsic absorption in excess to the 
Galactic value. 

• We detect (at more than 90% confidence level) the 
Fe Ka line in three objects. One of them is well de- 
scribed by a pure reflection model plus a Gaussian 
line at 6.4 keV rest-frame. This, the large equiva- 
lent width of the Fc line (although with large uncer- 
tainties) and diagnostics based on lines ratios from 
the optical spectrum support the hypothesis that 
this particular source is a Compton thick AGN. 

• We find four radio-quiet Type-2 QSOs. Their spec- 
tral energy distribution is well reproduced with a 
Seyfert-2 composite spectrum. 

• We confirm that in order to have a less biased sam- 
ple of AGN it is crucial to complement the opti- 
cal spectral properties with the X-ray informations 
(Lx and Nh), since many apparently normal galax- 
ies in the optical band are instead absorbed AGN. 

This is the first work on the X-ray spectral properties 
of the AGN in the COSMOS survey. We remark that 
once the XMM-Newton observations will be completed 
and the planned spectroscopic follow-up finished, we will 
be able to analyze the X-ray spectral properties of the 
AGN on a much larger sample and compare them with 
the properties of the AGN/host-galaxies at almost all the 
wavelengths. 
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TABLE 5 

Comparison between optical and X-ray classifications 
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